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ABSTRACT 

HST ACS images reveal blue cores in four E+A, or post-starburst, galaxies. Follow-up spectroscopy 
shows that these cores have LINER spectra. The existence of LINERs, consistent with those in many 
elliptical galaxies, is yet one more piece of evidence that these post-merger, post-starburst, bulge- 
dominated galaxies will evolve into normal ellipticals. More interestingly, if LINERs are powered by 
low-luminosity AGN, their presence in these E-f As suggests that any rapid growth phase of the central 
black hole ended in rough concert with the cessation of star formation. This result emphasizes the 
importance of E-t-As for exploring how the evolution of black holes and AGN may be tied to that of 
galactic bulges. 

Subject headings: galaxies: evolution — galaxies: active — galaxies: starburst 



1. INTRODUCTION 

The strong correlation between black hole mass 
and galaxy bulge velocity dispersion (Af, — as] 
|R:rraresc & Merritt 2000; Gebhardt et al. 2000} suggests 
one of two things. Either early-type galaxies do not arise 
from dissipative mergers or there is a connection — per- 
haps causal — between the small-scale physics of black 
hole (BH) growth and the large-scale physics that orga- 
nizes the host galaxy morphology, kinematics, and stellar 
populations during and after the merger. Given the ob- 
servational evidence that gas-rich mergers do occur and 
that at least some produce pressure-supported, bulge- 
dominated remnants, it remains to understand how the 
processes that drive the evolution of the smallest and 
largest galactic scales are related. 

A key to resolving this question is to identify galax- 
ies undergoing large-scale transitions via mergers and to 
consider the properties of their cores. On-going merg- 
ers are too complicated to provide clear answers, while 
their likely remnants are too far removed from the merger 
event. E-l-A galaxies, which have post-starburst spectra, 
frequent tidal features, and the kinema tic and morpho- 
logical signatures of early type galaxies fZabludoff^t_^ 
R996; Norton et al. 2001; Chang et al. 2001; Yan g et al. 
l2004f) . are true transitional objects and thus plausible 
test cases. In this Letter, we explore whether there is 
evidence for a central BH/AGN in nearby E-t-As, and, 
if so, whether the core is consistent with those of early- 
type galaxies and evolving in concert with the galaxy as 
a whole. 

2. DATA: IMAGING AND SPECTROSCOPY 

Our sample is a subset of the 20 nearby (0.06 < z < 
0.12) E+As (defined to have strong Balmer absorption 
Hnes, (H/37(5) > 5.5 A, but httle or no [O II] emis- 
sion, EW[0 III < 2.5 A) in the L as Campanas Red- 
shift Survey ( Zabludoff 'etal] I1996D . Initial inspection 
of our HST WFPC2 and ACS B and R images^ reveals 
at leas t six galaxies with bright, blue, almost stellar-like 
cores (|Yang et al.ll200(^ . In Fig. we show the color 
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profiles for three of these galaxies (EA06, EA16, EA17) 
and for another blue-core galaxy (EAOIB) that was ob- 
served serendipitously. EAOIB is the dis turbed compan- 
ion g alaxy of EAOlA (originally EAl in lZabludoff et al] 
11996ft that we have now determined spectroscopically 
to lie at the same redshift as EAOlA and also to have 
an E-l-A spectrum. The EAOlA-B system is the first 
known binary E-l-A system and provides additional ev- 
idence that the E4-A phase of galaxy evolution is trig- 
gered by galaxy-galaxy interactions. The central blue 
core of each galaxy has A(i? — i?) ~ 0.3 with respect 
to the outer galaxy and a characteristic size of 1 kpc. 
With the exception of the blue cores, the morphologies 
of these galaxies are consistent with elliptical/SO galaxies 
(EAOIB, EA06, EA16) or early spiral galaxies (EA17). 

The blue cores, while too faint to extract independently 
from models of the underlying galaxy light, raise the p os- 
sibility of a central AGN fe.g.. lO'Connell et alJl2005l) or 
of a luminous, compact star cluster fe.g.. ICoHna et alJ 
I2002D . We therefore obtained long-slit spectra of the four 
blue-core E-|-As (and EAOlA) with the Inamori Magellan 
Areal Camera and Spectrograph (IMACS) on the Mag- 
ellan Baade 6.5 meter telescope between April 3 and 5, 
2005. We used the 300 line mm~^ grism and 0.7" wide 
slit, which resulted in a dispersion of 1.34 A pixel" ^ and a 
spectral resolution of 4.2 A over the wavelength range 
4000-8000 A. We took three 600s or 900s exposures for 
each galaxy to remove the cosmic rays. The frames were 
processed and flux-calibrated in the standard manner 
with IRAF. The following analysis uses one-dimensional 
spectra extracted from a 0.7" x 2.0" aperture. No X-ray 
data are yet available for these sources to help discrimi- 
nate between the possible origins of the blue cores. 

3. NEBULAR DIAGNOSTICS OF NUCLEAR ACTIVITY 

Optical nebular emission lines provide key diagnostics 
of the state of the interstellar medium and the source 
of ionization (AGN or hot stars). Measuring these lines 
presents a particular challenge in E-|-As because they are 
weak (by definition) and superposed on a complex stellar 

3 (B, R) corresponds to {F435W, F625W) for the ACS imaging 
(EA06, EA16, and EA17) and to (F439W, F702W) for the WFPC2 
imaging (EAOl), respectively. 
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Fig. 1. — {B — R) color profiles of four E+A galaxies with blue 
cores. We note the bulge fractions (bulge-to-total light) and Sersic 
indices n (from a si ngle Sersic r^/" p rofile fit) derived from the 
GALFIT algorithm IPeng et aT] 120021) . The spheroid-dominated 
morphologies and post-starburst spectra are consistent with these 
galaxies being in transition from late- to early-types. 
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Fig. 2. — Examples of our continuum and line fits. The leftmost 
panel shows the data (dark) and the best fit continuum model 
(light, offset for clarity). The two panels on the right zoom in on 
the H/3 and Ho lines. The upper lines show the data (dark) and 
the continuum model (light). The lower lines show the continuum 
subtracted spectrum (dark) and our fit to the nebular lines (light). 
All three panels have the same vertical scaling. 



continuum. Conventional continuum subtraction tech- 
niques are inadequate, especially for the Balmer emission 
lines, which can be swallowed by 5 A or more of stellar 
absorption. Instead, we fit a stellar popul ation model 
to the continuum following iTremonti et alJ (|2006). We 
use a library of template spectra drawn from the stellar 
population models of ^Bruzual & Chariot ( 2003) , which 
have a spectral resolution of ~ 3 A. We use templates 
spanning a range of ages (1, 10, 25, 50, 100, 250, 500, 
750 Myr, 1, 1.5, 3, 6, 13 Gyr) and metallicities (0.4 Zq, 
1 Zq, 2.5 Zq). To construct the best fitting model, we 
perform a non-negative least squares fit with dust atten- 
uation modeled as an additional free parameter. 

After subtracting the best-fitting stellar population 
model of the continuum, we fit the nebular emission 
lines. Because we are interested in recovering very weak 
nebular features, we adopt a special strategy: we fit all 
the emission lines with Gaussians simultaneously, requir- 
ing that the Balmer lines (H7, H/3, Ha) have the same 
line width and velocity offset, and likewise for the for- 
bidden lines ([O III] AA4959,5007, [N II] AA6548,6584, 
[S II] AA6717, 6731). This procedure minimizes the num- 
ber of free parameters and allows the stronger lines to 
constrain the weaker ones. The measured Ha emission 
line equivalent widths range from 1 to 5 A for Ha, and 
from 0.2 to 1.2 A for H/3. Measurements of lines this 
weak are only possible because of the high S/N of the 
spectra (10 - 40 per pixel) and our continuum subtrac- 
tion techniques'*. We show examples of our continuum 
and line fits in Figure El 

To obtain error estimates for the line fluxes that in- 
clude errors in the continuum subtraction, we adopt a 
bootstrap approach. For each spectrum, we take the dif- 
ference between the continuum and the data as a measure 
of the error. We randomly re-sample the errors in bins of 

The ITremonti et ahj 12009) code allows the line amplitudes 
to be negative or positive, hence there is no overall bias towards 
detecting lines in emission. 



500 A, add them to the model of the continuum and neb- 
ular lines, and fit the resultant spectrum with our code. 
We do this 1000 times per galaxy, and use the spread in 
the measured emission lines as our estimate of the error. 
The technique does not address possible systematic er- 
rors in the stellar population models. Therefore, we per- 
form an identical fitting procedu re, substi t uting the syn- 
thetic spectra of Gonzalez-Dclg ado et al.l ll2005|) for the 
empirical ones of Bruzual fc CharlotI 1)2003 ^ . The result- 
ing differences in line fluxes are small in EAOIA, EAOIB, 
and EA17, but exceed the random errors in EA06 and 
EA16 (see Fig.|3). 

The weak nebular emission lines extracted from our 
spectra provide a key AGN diagnostic heretofore unavail- 
able foi_jnos^^+^As^^bllowing Baldwin et al. (198^ 
and iVeilleux fc Osterbroclj <)1987j) . we plot the flux ra- 
tios, [O III] A5007/H/3 versus [N II] A6584/Ha, in Fig- 
ure 13 (widely referred to as the BPT diagram). For com- 
parison, we plot 254,548 emission line galaxies from the 
Sloan Digital Sky Survey (SDSS) Data Release 4. The 
SDSS data form two distinct and remarkably narrow se- 
quences in this plot, the first corresponding to gas ion- 
ized by massive, main-sequence stars (the star forming 
sequence), the second to gas ionized by other means (the 
Seyfert/LINER sequence). Kewlev et al. (2001) theoret- 
ically calibrate a limit (dashed curve in Fig. ^ above 
which galaxies cannot be explained by any possible com- 
bination of parameters i n a standard star forming model. 
iKauffmann et al.l 1)20031) use SDSS data to empirically 
calibrate a limit that more closely adheres to the star 
forming galaxy locus (solid curve in Fig. 

Our E-|-As occupy an interesting part of the BPT plot. 
EAOIA falls squarely on the locus of star forming galax- 
ies, suggesting that there is minor star formation that 
was undetected previously due to dilution by light from 
larger radii. In contrast, the other four sources lie well 
above th e star forming locus, eve n considering the uncer- 
tainties. IKauffmann et all 1)20031) quantify the position of 
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AGN in the BPT plot using a polar scheme centered at 
the point where the AGN leave the star forming galaxy 
locus. Galaxies are characterized by their distance D 
from the origin and by an angle $, which is zero along 
the positive [O III]/H/3 axis and increases clockwise. Our 
four blue-core E+As hav e D = 1 - 1.1 and $ = 30°-46°. 
iKaufFmann et all 1)2003(1 classify such galaxies as LIN- 
ERs. 

The nature of LINERs l|Heckmanll98(]D remains a puz- 
zle. LINERs may be low-luminosity AGN (LLAGN) aris- 
ing from low-rate or low radi ative efficiency accretion 
onto super-massive blac k holes ijHalpern fc Steineilll983l: 
iFerland fc Netzeilll983D or they may be the product of 
other mechanisms (see [FiliDPcnk^ 120031 and references 
therein) such as photo ionization b y young clusters during 
the Wolf-Rayet phase l|Barth k Sh ields 2000). For exam- 
ple, a 4 Myr old star cluster was identified as the domi- 
nant ionizi ng source of the LINER nucleus in NGC 4303 
ijColina et al. 2002), but subsequent Chandra imaging re- 
vealed a hard X-ra y point source indicative of an AGN 
l|Jimenez-Bail6n et a l. 2003). In the nuclei of our E-fAs, 
the Balmer absorption lines constrain the mass fraction 
of < 10 Myr old stars to be less than 0.1%, ruling out a 
major contribution from young stars. 

A more plausible alternative to ionization by 
an AGN is ionization by planetary nebula nuclei 
that appear in the late-phase evolution of interme- 
diate mass stars. (Indeed, some LINER galaxies 
have significant inte rr nediat e- age stellar populations 
llCid Fernandes et"an l2004t iGonzalez Dekado et all 
|2004|) .') Howeve r, PNe nuclei are ineffi cient producers 
of Ha photons l|Taniguchi et al.l 12000(1 . and the inter- 
mediate age stellar mass required to generate the Ha 
luminosities of our nuclear spectra (a few xlO^^ erg s~^) 
is large (~ 10^ Mq in the slit aperture). While such 
large masses are not ruled out, a LLAGN provides a 
more natural explanation. Indeed, many LINERs have 
nuclear activity difficult to explain via stellar processes: 
compact flat-spectrum radio cores or jets. X-ray point 
sources with power-law spectra, and broad Balmer 
emissi on |Ho 2004, and references therein). Maoz et al. 
(|2005f) found that 80% of LINERs with compact UV 
cores show significant variability, strongly linking 
LINERs to AGNs. Thus, the LINER spectra in our 
E-KAs are likely signposts of weak AGN activity, but we 
acknowledge that intermediate mass stars may provide 
an alternate or additional source of gas ionization. 

4. DISCUSSION 

Using spectral diagnostics, we discover LINERs in the 
centers of four E-l-A galaxies with blue cores. These 
LINERs have low luminosit ies similar to th ose in many 
nearby early- type galaxies l|Ho et al.1 11997(1 . The exis- 
tence of LINERs in E-f-As, which are typically post- 
merger, post-starburst, and bulge-dominated galaxies, 
is more evidence that E-|-As will evolve into early-type 
galaxies. 

If these LINERs are low-luminosity AGN, i.e., pow- 
ered by black holes, a more interesting aspect of our 
discovery is the apparently simultaneous end of star for- 
mation and significant BH growth. If we estimate the 
AG N bolometr ic luminosity from the [O III] luminos- 
ity ( Heckman et al. 2004) . and the BH mass and as- 
sociated Eddington luminosity from the M, — as re- 
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Fig. 3. — AGN diagnostic diagram. The grayscale shows the 
locus of SDSS galaxies. The dashed curve is the AGN division 
theoretically calibrated by Kewlev et al. (2001); t he solid c urve is 
the empirical division suggested by iKauffma nn ct al. (2003) . The 
diagonal dashed line denotes <I> = 25°, the fKauffmann et al. divi- 
sion between Seyfert 2s and LINERs. The points with error bars 
indicate our E+A sample, fit with the Bruzual & Chariot (2003) 
models. The square points show the effect of fitting the contin- 
uum with the models of Gonzalcz-Dclgado et al. 1 2005j) and are an 
indication of the systematic error in the measurements. 



lation l(Ferrarese fc Merrittll2000t IGebhardt erahl 12000^ 
and Norton et al.'s (2001) velocity dispersion, the black 
holes in these E-|-A's are currently radiating at only 0.1 - 
0.5 % of their theoretical maximum luminosity. Thus any 
period of strong AGN activity and significa nt BH growth, 
as mi ght occur during the gas-rich merger ((Sanders et alJ 
Il988tl . has now ended. The same is true for the star 
formation: after experiencing a burst within the last 
Gyr, these galaxies are now forming stars at rates of 
less than 0.01 — 0.08 Mq yr~^ (derived from the mea- 
sured Ha luminositi es assuming SFR (Mq yr~^) — L^ia/ 
1.27 X lO^^ergs s"^: lKennicuttill998ll . 

These results have implications for the intersection of 
the small-scale physics associated with BH evolution and 
the large-scale physics that determines the properties of 
the hos t galaxy's bulge. The coherence of the M,—(jb re- 
lation llFerra,rese Merrittjl200(l IGebhardt et a,ni200(Tl 
requires that these small- and large-scale processes are 
connected, either causally or through a third process, 
particularly in systems where M, or ub is changing 
by large factors in a short time. If these LINERs are 
LLAGNs that represent the end of BH growth, and the 
E+A phase marks the point at which the galactic bulge, 
enhanced via the galaxy-galaxy merger and resulting 
starburst, stops growing, the study of E+As will provide 
unique insights into the origin of this relation. 

If the cores of these E-|-As do harb or LLAGN, are the 
AGN really fading? .K auffmann et a l. (2003) find that 
powerful AGNs (Seyferts with large D and $ < 25°) lie 
in massive early-type galaxies with young stellar popula- 
tions, whereas weak AGNs (pure LINERs with large D 
and $ > 25°) inhabit normal early- type galaxies. Given 
that the blue core E-f As are evolving from late- to early- 
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type galaxies through the post-starburst phase, it is quite 
possible that their AGNs are also in transition (or fading) 
from the strong (Seyfert) to the weak (LINER) AGNs. 
Finding Seyfert nuclei among the youngest E+As would 
confirm this conjecture, but our sample is presently too 
small for this type of search. 

Is the apparent relationship between the small-scale 
(LLAGN-sized) and large-scale (bulge-sized) physics in 
E-t-A galaxies causal? Recent numerical simulations in- 
cluding BHs suggest that collisions between galaxies trig- 
ger an inflow of gas that causes a strong circumnuclear 
starburst and fuels BH accretion. A powerful quasar out- 
flow removes the gas from the inner region of the merger 
remnant, quenching the star formation on a re latively 
short timescale (~1 Gyr;lSprin gel et al.'''2005a'b'). These 
models ultimately produce elliptical galaxies within a few 
Gyr and, unlike past simulations without AGN- feedback 
ijMihos fc Hernauisti IT991 11996(1 . effectively predict a 
true E-l-A phase in which the merger remnant has no lin- 
gering star formation. It is possible, though not proved, 
that the LINERs in our blue core E-t-As are the relic 
(fading) AGNs that once caused the truncation of star 
formation in these post-starburst galaxies. If so, their 
presence suggests a path by which early-type galaxies 
arrive on the Af, — as relation and also solves the long- 
standing puzzle of what mechanism during the merger 
ends the star formation in E-|-As. 

What is the occurrence of LINER spectra in the gen- 
eral post-starbu rst galaxy population? In recent work, 
lYan et all l)2005D investigate the emission line properties 
of a large sample of galaxies drawn from the SDSS. They 
find that AGN are characterized by large [Oil] /Ha EW 
ratios and conclude that post-starburst samples defined 
with an [Oil] EW cut (such as ours) will be biased against 
AGN. However, our high S/N spectroscopy reveals that 
even E-|-As with negligible emissio n (1-5 A a t Ha) have 
line ratios characteristic of AGN. lYan et all l(200!Tf l are 
able to use the BPT diagram (i.e., our Figure 3) to clas- 
sify only 40% of their post-starburst sample owing to the 
weakness of the emission lines. Of the galaxies they can 
classify, ~ 90% harbor AGN. They suggest that many of 
the weaker-lined objects may have LINER-like spectra, 
which is consistent with our results. 

A final puzzle is the nature of the blue cores, which first 
led us to investigate these galaxies. While our ground- 
based spectral resolution is not ideal for addressing this 
question, the low emission line luminosity of our LINERs 
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and the lack of broad emission lines suggest that the blue 
light is not from an AGN continuum. It may instead arise 
from a (fading) circumnuclear starburst, which grew in 
concert with any nuclear activity. Blue cores are common 
in early-type galaxies at higher redshift (z > 0.5), when 
field spheroids were assembling. For example, 30% of the 
morphologically-selected elliptical galaxies in the Hub- 
ble Deep Field North have color inhomogeneities, mostly 
due to blue cores (Mcnantcau et al. 2001). Trcu et aQ 
(|2005(1 find that - 15% (2/14) of blue-core spheroids 
in the Great Observatories Origins Deep Survey North 
(GOODS-N) field have X-ray detections {Lx > 10^^ ergs 
s~^), suggesting the presence of AGN. Our blue-core, 
LINER E-|-As may be local examples of a phenomenon 
common at high redshift. 

5. CONCLUSIONS 

We identify four E-fA galaxies with blue cores, which 
are revealed by our follow-up spectroscopy to have 
LINER spectra. The existence of LINERs, similar 
to those in elliptical galaxies, is more evidence that 
E+A galaxies, with their post-starburst spectra, post- 
merger, gas-poor, bulge-dominated morphologies, and 
pressure-supported kinematics (e.g. , .No rton et al. 200^ 
iGhang et atri200H lYang et all HSjT evolve into nor- 
mal early-types. More interestingly, if LINERs are low- 
luminosity AGN, their presence in E-|-As suggests that 
any rapid growth phase of the central AGN has ended in 
rough concert with the star formation and therefore that 
the evolution of the black hole is tied to that of the galac- 
tic bulge. What is not clear from our work is whether 
the coupling between AGN and bulge evolution is causal, 
as is suggested by some theoretical models incorporat- 
ing AGN feedback (Soringel et al. 2005a b). The study 
of a large sample of E-|-As, including an investigation 
of any correlation between AGN strength and the time 
elapsed since the starburst, could provide a test of the 
AGN-feedback hypothesis, as those models predict that 
the black hole accretion rate peaks shortly after the star- 
burst and declines quickly as the merger remnant ages. 
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